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We have performed a systematical calculation of the pion (n, K, rf) and heavy pseudoscalar meson 
S-wave scattering lengths up to 0(p A ) in the chiral perturbation theory in the heavy quark symmetry 
limit. With the three scattering lengths from the lattice simulations as input we estimate the 
unknown low-energy constants. Then we predict all the other unmeasured scattering lengths. The 
analytical expressions and predictions may be helpful to future investigations. Especially we note 
that the DK scattering length is positive. Therefore their interaction is attractive, which helps to 
lower the mass of the "bare" charm-strange scalar state in the quark model through the couple- 
channel effect. 
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I. INTRODUCTION 

The discovery of the narrow D s j(2317) has inspired 
heated discussions of its structure in the past six years 
[E 0, HI ■ The possible interpretations include the chiral 
partners of D s [A], P-wave excited states of D S (D* S ) 
couple-channel effects between the cs state and DK con- 
tinuum [fl conventional cs states @, 11] , four-quark states 
0, [T3JlIl [H, [H, [H, [ill etc. (for a detailed review see 
RefTlll)- 

Considering the large contribution of the S-wave DK 
continuum, the mass of D*j(2317) agrees well with the 
experimental value, which is indicated by the coupled- 
channel effect @ and the QCD sum rule (QSR) approach 
[l~7l ] . Therefore, D*j(2317) is very probably a conven- 
tional cs state with J p — + . Since D s j(2317) strongly 
couples to the DK channel, the study of the DK inter- 
action is very interesting. 

The scattering length is an important observable, 
which encodes the information of the underlying interac- 
tion. For example, a positive scattering length suggests 
there exists attraction in this channel. In this work we 
investigate the scattering lengths in the pion-heavy me- 
son channels. Here the pion denotes tt, K or rj while the 
heavy mesons are the pseudoscalar charmed or bottom 
mesons. 

There are several works on these scattering lengths in 
the literature. The S-wave DK scattering length was pre- 
dicted to be 5 ± 1 GeV -1 in a unitarization model llal. 
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a D^ = 0.29(4)m~ 1 were extracted from a lattice cal- 
culations of the scalar form factors in the semilcptonic 
decays 19]. Recently, a lattice study on the scattering 
lengths of the light hadrons with the charmed mesons 
and charmonia were performed in full QCD (20| . Un- 
fortunately their pion masses are still quite large. The 
scattering lengths of heavy mesons and Goldstone bosons 
were discussed to the next leading order in chiral per- 
turbation theory recently [2lj . In this paper, we will 
go to the next-next-leading order of the chiral expansion 
and consider the important loop corrections to the elas- 
tic pion-heavy meson scattering at threshold in the heavy 
quark limit. 

It is known that the chiral perturbation theory works 
well for the light pseudoscalar meson systems. For 
the meson-baryon interactions, the heavy baryon chi- 
ral perturbation theory (HByPT) was proposed so that 
a systematic power counting rule exists. The scat- 
tering lengths in ttN, KN and other channels have 
been investigated to high orders within this framework 
[22I I23I l24l [25[ | . It was observed that the chiral expansion 
in the SU(3) case converges well only in few channels. 

For the heavy mesons, heavy quark symmetry [26] im- 
poses that the ground states D, D s , D* and D* belong 
to the same doublet. Similar to the meson-baryon case, 
the heavy meson chiral perturbation theory (HMyPT) is 
a useful tool to study their interactions [27], . In the 
heavy quark symmetry limit, the recoil order corrections 
are neglected. One may use a power counting rule sim- 
ilar to that in HBxPT to include the chiral corrections 
order by order. We will adopt the HMxPT formalism to 
calculate the scatting lengths up to 0(p 3 ). Such a study 
may be helpful to test the convergence of HMxPT in the 
scattering processes. 

In our formalism, the S-wave scattering length is de- 
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fined through 

vn 

T th = 8tt(1 + — )o (1) 

where is the threshold T-matrix element and m (M) 
denotes the light (heavy) meson mass. 

This paper is organized as follows. In Sec. UH we 
present the chiral Lagrangians for the calculation at 
threshold. In Sec. 1111} we give the threshold T-matrices 
for the S-wave elastic scattering of the Goldstone bosons 
and charmed mesons up to the third order. According to 
the heavy quark symmetry, these T-matrices are also ap- 
plicable to the pion-bottom meson case. In Sec. IIV1 we 
estimate the values of the scattering lengths. The final 
section is the discussion. 



II. LAGRANGIANS 



flavor indices is implicit. The heavy fields have been re- 
defined with a factor \[M and have the mass dimension 
3/2. 

For the calculation at threshold, we need the following 
0(p 2 ) Lagrangian 

C%\ = co(HH)ti( X+ ) - c^Hx+H) 

-c 2 {HH)tr(u u ) - c 3 {Hu u a H). (9) 

The signs of these terms are consistent with those in Ref. 
& 

Similar to the chiral Lagrangians for the meson-baryon 
systems (30L f3l| , there are many terms in the third order 
Lagrangian. For simplicity, we write down only the piece 
relevant to our threshold calculation: 



The lowest chiral Lagrangian for the light pseudoscalar 
mesons reads 



(2) 



chiral limit, is the axial vector field with the definition 



where / « 92.4 MeV is the pion decay constant in the 
the axial vector field with the 

;U f A£}, £ = exp(z0/2/), (3) 
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and 



X± = &X&±ZxZ, X = diag(m^, ml, 2m 2 K -ml). (5) 
The doublet of ground state heavy mesons reads 

5 = 7 o F t 7 o = (p*V + iP t 75 )iM ) 
P=(D°, D+,D+), P* = (D°*,D+*,Df*)^ (6) 

or 

P={B-,B°,B° S ), P* = (B-*,B°*,B°*)^, (7) 

where = (1,0,0,0) is the heavy meson velocity. The 
leading order chiral Lagrangian for the heavy mesons in 
the heavy quark symmetry limit is 

£hI = -((id H)H) + (HTqH) + g{Hu^ l5 H), (8) 

where T M = , d^] is the chiral connection, (...) means 
the trace for Gamma matrices and the summation over 



C% = k(H\x-,uo]H). 



(10) 



III. THRESHOLD T-MATRICES 



We omit details and present the explicit expressions 
for threshold T-matrices to 0(p 3 ) of the chiral expansion. 
The numerical results will be estimated in the next sec- 
tion. At the leading order, we have Weinberg- Tomozawa 
terms 
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with K = (K+,K a ) T , K = {K°,K-) T . The super- 
scripts of T-matrices denote the total isospin. Here we 
have replaced / with the renormalized decay constants 
/tti Ik and f v . The resulting 0(p 3 ) corrections will be 
taken into account later. These expressions are consis- 
tent with the results in Ref. [2l[ considering the different 
overall sign for T. 

It is convenient to define two combinations of the low 
energy constants (LECs) at the next-leading order 



Ci = 8co-4 Cl + 2c 2 + C3, 
C = 8c + 4ci + 2c 2 -c 3 . 
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(2) 

From C H( p, we get the following T- matrices 



T, 



L<h ~ 2/ 

(1) 



(o) 

DK 



-Ci, 



DK Ik DK 7|~ 



Td s k 

(3/2) 



-Ci, 



T, 



f2 ^l, ^D S K — f-! 

l Dtt — fl °1) J Dtt — /2 L/ l) 
m 2 

2r> s 7T = jjk(Ci + Co) 
T Dv = 3)7 [(2d + C )to 2 + 4ci ( TO 2 _ m 2 J] ; 

= [(7d - C )to 2 - 16 Ci (to 2 _ m 2)] . (13 ) 

The effects due to the decay constant renormalization are 
beyond the order we are considering. 

At the third order, we have the non-vanishing correc- 
tions from the loop diagrams in Fig. [TJ 
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In deriving the above analytical expressions, we have 
used the dimensional regularization and minimal subtrac- 
tion scheme. Here we retain the chiral logarithm terms. 
A ~ 1 GeV denotes the scale of chiral symmetry break- 
ing. In these T-matrices, the corrections proportional to 
g 2 come from intermediate vector meson contributions 
(the last six diagrams depicted in Fig[T]) and the diver- 
gences from these diagrams cancel out. The fourth dia- 
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gram generates imaginary parts for T D X K , T D K , Td v 
Td s 7j- 



and 



FIG. 1: Non- vanishing loop diagrams in the calculation of meson-heavy meson scattering lengths to the third chiral order in 
HM/yPT. Dashed lines represent Goldstone bosons while solid (double) lines represent pseudoscalar (vector) heavy mesons. 
The fourth diagram generates imaginary parts for kaon-heavy meson and eta-heavy meson scattering lengths. 



From Cjj,, we obtain the counter terms: and the SU(3) relations (m„ = = ni s ) 



T& = T D . T =$[TV+Tgl], (29) 
T« =0 , 2^ = 16^, T DsK =T^=T^l (30) 

T _ lr T (3/2) 9T (l/2)i _ lr T (l) , T (0) i /oi\ 

T Ds7r + T DbV = \ [2T§1 2) + Tgf 3 ] + T Dri . (32) 
These relations are very useful to cross-check the results. 

( 25 ) IV. NUMERICAL RESULTS 
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where n r — n r (X) = n 



|-p is the renormalized LEC 
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One may verify that the above threshold T-matrices 
satisfy the crossing symmetry: 



There are four parameters to be determined in our T- 
matrices: c\, C%, Co and k t . We extract ci from the mass 
splitting between heavy mesons within the same doublet: 
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(28) 



where M is the heavy meson mass in the heavy quark 

o 

limit and chiral limit. We take the value M~ 1918 MeV 
for the charmed system (29|. Using M D = 1867.2 MeV, 
M Ds = 1968.5 MeV, M D . = 2008.6 MeV, and M D , = 
2112.3 MeV (H, we get a = 0.47 GeV" 1 . 
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Up to now, there does not exist any available experi- 
mental measurement of these scattering lengths. A recent 

(3/2) (1) 

lattice simulation yielded a D ' v = —0.16(4) fm, a D g = 
-0.23(4) fm, a Dsn = 0.00(1) fm and a DsK = -0.31(2) 
fm [2fj|. In order to determine the other three LECs, we 
use these values as inputs. 

From a^ 2 \ ar> s if and a_D s u-, we get G\ = —0.63 
GeV-\ C = 5.64 GeV^ 1 and n r = -0.33 GeV" 2 . 
We have used = 139.6 MeV, m K = 493.7 MeV, 
m v = 547.8 MeV, f n = 92 A MeV, f K = 113 MeV, 
f v = 1.2/tt and A = 47r/ T . The coupling constant 
q = 0.59 ± 0.07 ± 0.01 was determined from the width of 

D*+ mm- 

We present our numerical estimations for T-matrices 
and scattering lengths in Table [U For comparison, wc 
also list several scattering lengths appearing in the liter- 
ature. 

The heavy quark symmetry works better for bottom 
systems. We would like to estimate the S-wave scattering 
lengths for the elastic pion-bottom meson interactions as 
well. 

In the bottom case, one extracts c\ — 0.39 GeV" 1 
with the mass 32] M B = 5279.3 MeV, M B , = 5366.3 
MeV, M B * = 5325.1 MeV, M B * = 5412.8 MeV and 

M= \{M B + M Bs ) = 5322.8 MeV. The heavy quark fla- 
vor symmetry tells us that we may still use the above 
determined values for C\, Cq and k t . The coupling 
constant g from a recent unqucnched lattice result is 
0.516(5) stat (33) c fci rai (28)pert(28) disc [35). We use g = 
0.52 to estimate the values. The calculated scattering 
lengths are presented in Table ITT1 

In principle, the scale dependence from the chiral loop 
is compensated by LECs. By varying the scale A = 4^/^ 
to A = m p = 775.49 MeV, one may check the ana- 
lytic independence on it. By repeating the above pro- 
cedure, we get d = -0.75 GeV" 1 , C = 5.76 GeV" 1 
and k t = —0.21 GeV -2 at this lower scale. We show the 
numerical estimations for the charmed systems in Table 
IIIII The scale independence is illustrated if one compares 
the results with those in Table U 



V. DISCUSSIONS 

We discuss the results at the scale A = Airf^. By com- 
paring the values in Table |TT] and [H one notices the dif- 
ference between the results of pions-bottom mesons and 
those of pions-charmed mesons is small. Larger differ- 
ence will appear after including the recoil corrections. 
We focus mainly on charmed systems in the following 
discussions. 

The positive sign of a^^, a^ K , , a^ 2 -* and a Bn 
indicates that the interactions for these channels arc all 
attractive. For the isoscalar D — K channel, the attrac- 
tion is relatively strong. However, further exploration of 
the phase shifts of the elastic DK scattering is required 
in order to answer whether the DK interaction is strong 



enough to form a bound state such as a DK molecular 
state. 

The result for the isoscalar D — K channel is very in- 
teresting. Since no quark pair annihilation occurs in this 
channel, the relatively large scattering length implies that 
the possibility to form a four-quark resonance state is not 
excluded. If the strong attraction is confirmed in future 
studies, this channel is certainly worthwhile further ex- 
ploration. 

The D — 7/ channel is also interesting. Twenty years 
ago, the formation of 77-mesic nucleus [36j was proposed 
based on the observation that rj is neutral and the strong 
force between rj and N is attractive. Later the possibility 
of the 77-hypernuclei was also discussed 0- I* 1 ^ e ^- 
the bound state of a D meson and 208 Pb was studied. 
Now, with the observation that both Drj and 77^ are at- 
tractive, one may guess the formation of a 77 — I?— heavy 
nucleus state is also possible. However, such a bound 
state would be very difficult to be detected experimen- 
tally Hi- 

Contrary to the D — ij channel, the scattering length 
aD gTI is negative, which indicates the interaction is repul- 
sive. In order to get a positive cld^ and nearly vanishing 
o-D,-n ~ 0, one requires C\ > 1.4 GeV -1 and Co < 3.6 
GeV" 1 . 

From Table [H we note that the chiral expansion con- 
verges well in four channels: isoscalar DK, isovector DK, 
D S K and D S K. The convergence of T^ K , Tur, and Izj s?) 
will be manifest when the 0(p 4 ) corrections are included. 
Although the isoscalar D — K channel is very interesting, 
there is a critical problem of convergence. The contribu- 
tion from the second order is very large. Future accurate 
determination of LECs is probably helpful to diminish 
this contribution. 

Future improvements may be made for the present cal- 
culation. We have used the recent lattice results in full 
QCD to estimate the LECs where the lattice pion mass 
is large [2(J. Due to the lack of data, the deduced LECs 
have large uncertainties. Since no J p = + resonance 
close to the D-k threshold has been reported, this chan- 
nel is ideal to determine the LECs. To test HMxPT 
in the scattering problem, the experimental determina- 
tion of pion-heavy meson scattering lengths is strongly 
expected. On the other hand, further lattice simulations 
can provide more reliable inputs. We hope our chiral cor- 
rections may be useful when performing chiral extrapo- 
lations. 

One may also improve the analytical expressions by 
including 1/M corrections and extending the calculation 
to higher orders in future exploration. In addition, the 
effects of the nearby + resonances may affect the predic- 
tions. One can consider them using the non-perturbative 
methods. 

Finally, we mention that the scattering lengths of pions 
with heavy anti-mesons. They are easy to get through 
the C parity transformation 

T (i) _ T (i) T (i) T (i) _ T (i) ( ~ A s 

HK HK ' HK ffif > 1 Hn/, 1 ~ 1 Htt/ v ^ 
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o( P 3 ) 
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Scattering lengths 


Other results 


DK 





9.4 


-1.5 + 5.6* 
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0.0 
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0.34 + 0.17; 




TD s ri 


0.0 


-13.8 


0.5 + 11. li 


-13.3 + 11. li 


-0.41 + 0.35; 





TABLE I: Threshold T-matrices for the elastic scattering of pions and charmed pseudoscalar mesons in unit of fm with the 
scale A = 4irf x . For the results from a unitarized method, one may consult Ref. [2l| . 





0(p) 0( P 2 ) 


o(p 3 ) 
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Scattering lengths 
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T (l) 

BK 





9.4 


-1.4 + 5.6i 


8.0 + 5.6; 
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T (l) 
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-0.42 
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7.6 


21.2 
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1.16 
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-7.6 


-2.4 


-0.1 


-10.1 


-0.37 




Tb s k 
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-4.3 
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T (l/2) 
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6.4 


0.25 


0.37 + 0.37 [19] 







1.1 


-1.1 
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0.00 




Tb v 
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0.9 + 5.5; 


10.3 + 5.5i 


0.37 + 0.20; 




Tb,ti 





-12.8 


0.3 + 11. li 


-12.4 + 11. li 


-0.45 + 0.40i 





TABLE II: Threshold T-matrices for the elastic scattering of pions and bottom pseudoscalar mesons in unit of fm with the 
scale A = 4tt/ x . Here B = (B°, B-) T . 



where / is the total isospin and H (H) denotes the heavy 
meson (anti-meson). 
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-7.6 


-2.8 


-0.7 


-11.2 


-0.35 


T (0) 
DK 


7.6 


21.7 


2.1 


31.4 


0.99 


Td s k 


-7.6 


-2.8 


0.7 


-9.7 


-0.31 (input) 


D K 


7.6 


-2.8 


-1.9 + 8.3z 


2.9 + 8. 3i 


0.09 + 0.27i 


r (3/2) 


-3.2 


-0.3 


-0.8 


-4.3 


-0.16 (input) 


T (l/2) 
1 Dtt 


6.5 


-0.3 


0.2 


6.3 


0.23 


Td s tv 





1.1 


-1.1 


0.0 


0.00 (input) 


Tdt) 





9.7 


1.2 + 5. 5i 


10.8 + 5.5i 


0.33 + 0.17i 


Td s ti 





-14.5 


0.5 + 11. li 


-14.1 + 11. li 


-0.44 + 0.35i 



TABLE III: Threshold T-matrices for the elastic scattering of pions and charmed pseudoscalar mesons in unit of fm with a 
lower scale A = m p . 
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